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Abstract. G. Giuga conjectured that if an integer n satisfies k n ~ l = — 1 

k=\ 

mod n, then n must be a prime. We survey what is known about this 
interesting and now fairly old conjecture. 

Giuga proved that n is a counterexample to his conjecture if and only if each 
prime divisor p of n satisfies (p — 1) | (n/p — 1) and p | ( n/p — 1). Using 
this characterization, he proved computationally that any counterexample 
has at least 1,000 digits; equipped with more computing power, E. Bedocchi 
later raised this bound to 1,700 digits. By improving on their method, we 
determine that any counterexample has at least 13,800 digits. 

We also give some new results on the second of the above conditions. This 
leads, in our opinion, to some interesting questions about what we call Giuga 
numbers and Giuga sequences. 
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1 Introduction 


In 1950, G. Giuga formulated the following conjecture ([3]). 

Conjecture. For each integer n it is true that 

n— 1 

n is prime s n := ^ k n ~ l = —1 rnodn. 

k=\ 

Fermat’s little theorem says that if p is a prime, then k p_1 = 1 nrodp for 
k = 1,... ,p — 1. Therefore, for each prime p, s p = — 1 rnodp. The question 
becomes: 

Does there exist a non-prime n such that s n = — 1 rnodn? 

This question has resisted solution for more than forty years. After surveying 
what is known about the conjecture, we will give several new results here 
which might suggest directions of further investigations. 

The key to dealing with Giuga’s conjecture is the following theorem, which 
was proved by Giuga in his original paper. A proof can also be found in [5]. 
For the sake of completeness, we give the proof here. 


Theorem 1. s n = — 1 rnodn if and only if for each prime divisor p of n 
we have (p — 1)1 { n /p — 1) and p | ( n/p — 1). 


Proof. It is well-known and an easy consequence of considering residue 
classes (see [5], p. 16) that for a prime p, we have 

h n -1 = f - 1 mod p if (p - 1) | (n — 1), 

^ ~ \ 0 mod p if (p — 1) / (n — 1). 

Therefore, for each prime divisor p of n with n = p ■ q, we get 


n—1 

fc=i 




k =1 


—q mod p if (p — 1) | (n — 1), 
0 mod p if (p — 1) / (?r — 1). 


(1) 


Assume that s n 
we have 


= — 1 mod n. Then for each prime divisor p of n, n = p ■ q, 

_ 1= f-9 modp if (p — 1) | (n —1), ( ) 

~ \ 0 mod p if (p — 1) / (n — 1). 
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This is only possible if (p— 1) | (n —1) = q(p—l) + (q — l). So (p— 1) | (q — 1). 
It then also follows from (2) that — 1 = — q mod p, or p \ (q — 1). 

On the other hand, assume that p \ (q — 1) and (jp — 1) | (q — 1). It then 
follows from (1) that s n = — q mod p; since q = 1 mod p. we have that 
s n = — 1 mod p for each prime divisor p of n. Now, n must be squarefree: 
If it were not, then there would exist a prime divisor p of n with p \ q; this 
contradicts p \ (q — 1). Since each of the distinct prime divisors of n divides 
s n + 1, this is also true for n. In other words, s n = —1 mod n. © 

As noted in the proof, every counterexample to Giuga’s conjecture must 
be squarefree. Squarefree composite numbers which satisfy the first of 
these two conditions have been investigated in their own right: they are 
called Carmichael numbers. They were introduced by Carmichael in 1910. 
Carmichael numbers are of interest because they are “pseudo-prime” in the 
following sense (Korselt’s criterion, 1899): n divides a n — a for all integers a 
if and only if n is squarefree and p— 1 divides n/p— 1 for all prime divisors 
p of n. The Carmichael condition 

(p — 1) | (n/p — 1) for all prime divisors p of n 

is equivalent to the condition 

(p — 1) | (n — 1) for all prime divisors p of n. 

Note that any Carmichael number is odd. (Assume that n is even. It has 
at least one other prime factor p besides 2. Then the even number p — 1 
divides the odd number n — 1, which is a contradiction.) The smallest 
Carmichael number is 561 = 3 • 11 • 17. The next two Carmichael numbers 
are 1105 = 5 • 13 • 17 and 1729 = 7 • 13 • 19. It has only recently been proved 
that there are infinitely many Carmichael numbers (see [1]). 

In order to refer easily to the second condition as well, we will call any 
composite number n with p \ (n/p — 1) for all prime divisors p of n a 
Giuga number. As we saw in the proof of Theorem 1, any Giuga number 
is squarefree. Moreover, one can prove the following equivalence along the 

lines of that proof (with n— 1 replaced by p(n), where cp is the Euler (totient) 

ra-i 

function): n is a Giuga number if and only if k^ 71 ) = — 1 mod n. In his 

k= 1 

original paper, Giuga proved another equivalence: n is a Giuga number if 
and only if 

V--TT-e^V- 

i\ n p 


3 



This equivalence will be of great importance throughout this paper. We 
will give a slightly generalized version of Giuga’s proof below. The smallest 
Giuga number is 30: 1/2 + 1/3 + 1/5 — 1/30 = 1. The next two Giuga 
numbers are 858: 1/2 + 1/3 + 1/11 + 1/13 - 1/858 = 1, and 1722: 1/2 + 
1/3 + 1/7+ 1/41 — 1/1722 = 1. We do not know if there are infinitely many 
Giuga numbers. 

Giuga’s theorem can now be restated as 

Theorem. A composite integer n satisfies s n = —1 rnodn if and only if it 
is both a Carmichael number and a Giuga number. 

Giuga’s conjecture is that such a number cannot exist. 

The smallest odd Giuga number has at least 9 prime factors, since with a 
smaller number of prime factors the sum 1/pi + ... + 1 / p m — 1/n is smaller 
than 1. However, this lower bound for a counterexample increases dramati¬ 
cally if we take into account that it must also be a Carmichael number. Any 
Carmichael number n has the following property: If p is a prime factor of n, 
then for no k is kp + 1 a prime factor of n. (If it were, then we would have 
( kp + 1) — 1 = kp | (n — 1) and p | n, which is a contradiction.) So, for 
example, no Carmichael number has the prime factors 3 and 7 at the same 
time. This property was used by Giuga to prove computationally that each 
counterexample has at least 1000 digits. Later, E. Bedocchi ([2]) used the 
same method to prove that each counterexample has at least 1700 digits. 
We will describe the method in Section 2 of this article. We have been able 
to improve on this method by reducing the number of cases to be looked at 
and have shown computationally that any counterexample has no less than 
13,800 digits. 

We believe that an approach to prove or refute Giuga’s conjecture in general 
is to study Giuga numbers in more depth; we will do this in Section 3 of 
this paper. It turns out that there is much more structure to be studied if 
we drop the condition that the numbers p in the definition be prime. This 
leads us to the following definition. 

Definition. A finite, increasing sequence of integers, [ni,...,n m ], is called 
a Giuga sequence if 



The proof of the following equivalence is due to Giuga in the case that the ni 
are primes. 
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Theorem 2. A finite, increasing sequence [n±,... ,n m ] is a Giuga sequence 
if and only if it satisfies ni \ (n\ ■ ■ ■ rii-i ■ nj + 1 • • • n m — 1) for i = 1,..., m. 

Proof. Write n := n\ ■ ■ ■ n m and qi := n/n^. Note that the sequence is a 
Giuga sequence if and only if n | (q\ + ... + q m — 1). 

The “only if” part of the asserted equivalence now follows immediately from 
this. 

On the other hand, assume that w* | (q t — 1) or n| | (n — nf) for all i. 
Multiplying leads to n 2 | (n — n{) ■ ■ ■ (n — n m ). In evaluating this product 
we can drop all multiples of n 2 . We therefore get n 2 \ (n (q\ + ... + q m ) — 
n\ ■ ■ ■ nm) = n (qi + . ■ ■ + q m - 1). Therefore, n \ (q\ + ... + q m - 1), which 
means that [n \,..., n m \ is a Giuga sequence. © 

Note that each two distinct elements rq, nj in a Giuga sequence are relatively 
prime (n* | (n/rii — 1), but rij \ nfnf). When n is a Giuga number, it gives 
rise to a Giuga sequence (its prime factors), but in general it is conceivable 
that an integer n can have two different factorizations, both of which are 
Giuga sequences. However, we know of no example of this. We do know 
that there is an infinity of Giuga sequences; we will show this in Section 3 
of this article. 

In the same way, Carmichael numbers can be generalized to Carmichael 
sequences: A finite, increasing sequence, [a\,... ,a m \, is called a Carmichael 
sequence if [a l — 1) | (ai • • • a m — 1) for i = 1,..., m. Note that a Carmichael 
sequence has either exclusively odd or exclusively even elements and that 
its elements need not be relatively prime. As with Carmichael numbers, any 
two distinct factors ai,aj in a Carmichael sequence satisfy a* ^ 1 mod aj 
(or, equivalently, aj / (a* — 1)). Carmichael sequences occur in much greater 
profusion than Giuga sequences. At the end of Section 3 we will construct 
some infinite families of Carmichael sequences. 

Giuga’s conjecture would be proved if one were to show that no Giuga 
sequence can be a Carmichael sequence. This, in turn, would be proved if it 
can be shown that any Giuga sequence must contain two factors n^, nj with 
nj | (ni — 1). It might even be true that every Giuga sequence contains an 
even factor; Giuga’s conjecture would follow from this. We have found no 
Giuga sequence which consists of odd factors only, but this is probably a 
consequence of the size of the problem. 

In Section 4 of this article, we will give a list of open questions concerning 
Giuga sequences and Giuga’s conjecture. 
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2 Computing lower bounds for a counterexample 


As we have seen, any counterexample to Giuga’s conjecture must be a 
squarefree odd number with prime factorization n = q\ ■ ■ ■ q k such that 

(i) qi ^ 1 mod qj for all i.j, and 

(ii) 1/ q\ + ... + l/q k > 1. 

Giuga and Bedocchi used these two properties to compute lower bounds 
for a counterexample in the following way. For m £ JN, denote by p m the 
to- th odd prime. A finite set of odd primes, {q\, ..., q k }, is called normal if 
condition (i) obtains. For each to £ IN, let S m be the set of all normal sets 
with maximum element smaller than p m . For each S £ S m , S = {(/i,..., q k }, 
define the set 

F m (S) = {(j 1 !j • • • j <?fc, q k +h • • • > Qr} 

to be the smallest set of odd primes which contains S, and is such that 

r 

Qj > Pm and S U {qj} is normal for j > k, and J2 1 /Qj > 1- Let r m (S) be 

3 = 1 

the number of elements of T m {S). For example, 

T 6 ({3,5}) = {3,5, 17,23,29,47,53,... ,7919,7937}, r 6 ({3,5}) = 383 
and 

Ti({}) = {3,5, 7,11,13,17,19,23,29}, n({}) = 9. 

Now define the sequence (i m ) m6 jv by 

i m = rnin{r m (5) | S £ S m }. 

As Giuga observed, this sequence is non-decreasing; we will see shortly why 
this is the case. 

Now, the number of prime factors of a counterexample to Giuga’s conjecture 
exceeds i m for each to £ IN. Indeed, the prime factors form a normal set, 
and the subset S of the factors smaller than p m is a member of S m . Since any 
normal set of primes which contains S and satisfies condition (ii) above must 
contain at least r m (S) elements, we have that n has at least r m (S) > i m 
prime factors; this is true for any to. £ IN. So, any counterexample is 

im 

bigger than J} pj , and therefore has at least the same number of digits as 
3 =1 

this product. Giuga estimated ig > 361, this yields more than 1000 digits; 
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Bedocchi computed ig = 554, this yields more than 1700 digits. (Note that 
Giuga and Bedocchi used a slightly different definition for i m \ this is why 
their numbers differ from our numbers by 1.) 


To compute i m , one has to find r m (S) for all S £ S m . Since the number 
of elements of S m increases geometrically with m, the time needed to com¬ 
pute i m gets out of hand quickly. With our R4000 Challenge server and 
the symbolic manipulation package Maple, we were able, with considerable 
effort, to compute i\g = 825. At this point in time we started looking for 
a better algorithm, something which allows us to look at only some sets 
in S m , not all of them. Fortunately enough, we found just such an algo¬ 
rithm. It is based on the following observation. Consider a set S £ S m and 
the associated value r m {S). Now, S has at most two “successors” in the 
set namely S itself and the set S' = <SU{p m }. We will now show that 

r m+ i(S) > r m (S) and r m _|_i(S") > r m (S). In fact, there are two cases: 


Case (i): S U {p m } is normal. Then S has the two successors S and S' 
in S m . |_i. Also, we have p m £ T. m (S). However, p m 0 T. m+ i(S), but every 
other element of T. m (S) is contained in T m+ i(S). So, T m+ i(S) must contain 

at least one higher prime for the sum 1/(7 to exceed 1. Therefore, 

q£Tm + i(S) 

r m+ i(S) > r m (S). As regards S' , the set T. m {S) may contain primes which 
are congruent to 1 mod p m . These are missing in T m+ i(S'), since p m £ S'. 

For each of these we need at least one higher prime for the sum 1/q 

q£T m+1 {S ') 


to exceed 1. Again, > r m (S). 


Case (%%): S U {p m } is not normal. Then the only successor of S in S m+ 1 
is S itself. Also, T. m (S) = T m+ \(S)] the prime p m is not contained in either 
set. Therefore, = r m _(_i(5). 


This shows that the sequence i m is indeed non-decreasing. But it shows 
more: the values rk+i, fk+ 2 , ■ ■ ■ for all of the successors in Sk+i, Sk+ 2 , ■ ■ ■ 
of a given set S £ Sk do not fall below r^S). If we want to compute i m 
and already know an upper bound I > i m , then we do not have to look at 
any successor in the sets Sk+ 1 , • • •, S m of a set S £ Sk with rfc(S') > I. So 
the natural way to do this is to do it iteratively: Start with Ag := Si, and 
let Ak -|_i consist of the successors in Sk+i of all S £ Ak with rk{S) < I. 
Then i rn = rnin{r m (5) | S £ A m }. If I is close to i m , then this significantly 
reduces the number of sets to consider. 


The bound I can of course be chosen as the value r m (S) for some S £ S m . 
The iterative method saves the most time if one correctly guesses which 


7 



sets have low values. By looking at preliminary computational results, we 
discovered that the following rule seems to hold. Let L 5 := {5, 7}, and define 



Lk U {pk} 
Lk 


if Lk U {pk} is normal, 
otherwise. 


Then it seems that for m > 5, r m (L m ) = i m . We have no proof that this 
is always true, but having discovered that the sets L rn yield good upper 
bounds for i m , we employed our iterative method with these upper bounds 
to compute all values i m for m < 100 in Maple and later for rn < 135 in C. 
We always found that r m (L m ) = i m . 

It was, by the way, surprisingly difficult to translate what was a fairly 
straight-forward Maple program into C. While Maple handled the data 
structures we required (lists of sets of variable length) easily, it was a non¬ 
trivial problem to implement these in C. We gained a speed-up of a factor 
of around 5 (for m around 100), though. Even with this speed-up, the last 
case (rn = 135) took 303 cpu hours, and the “curse of exponentiality” makes 
further computation close to impracticable. We thank Gerald Kuch (now 
a graduate student at the University of Waterloo) for doing this conversion 
from Maple to C. 

Here are some of the i m (the first nine of these are also given by Bedocchi 
in [2]): i\ = 9, h = 27, *3 = 65, 14 = 114, = 127, = 202, i-j = 278, 
is = 323, ig = iio = in = 554, ii 2 = ii3 = hi = ii5 = ii 6 = 704, 
hi = ii 8 = 751, ii 9 = i 2 o = 825, ..., 149 = i .50 = 2121, ...,*74 = *75 = 2657, 
• • •, *99 = *100 = iioi = 3050, ..., *131 = *132 = *133 = *134 = ii 35 = 3459. 
iioo = 3050 implies that any counterexample to Giuga’s conjecture has at 
least 12055 digits, *135 = 3459 implies that any counterexample has at least 
13887 digits. 


As Bedocchi points out in [2], this method is inherently incapable of show¬ 
ing that Giuga’s conjecture holds for all integers: the set L 27692 is normal, 
has 8135 elements and satisfies J2 1/(7 > 1. Therefore, i m < 8135 for 

g£i27692 


all m > 27692. 



3 Giuga sequences 


Recall that a Giuga sequence is a finite, increasing sequence of integers, 
[77-1 ,..., n m ], such that 



When dealing with Giuga’s conjecture, we are mainly interested in Giuga 
sequences which consist exclusively of primes; we will call these proper Giuga 
sequences. However, Giuga sequences, proper or not, are interesting objects 
in their own right; in this section we will give some of their properties. 

We have computed all Giuga sequences up to length 7 and some of length 8. 
There is no Giuga sequence of length 2; one sequence of length 3 ([2,3,5]); 
two sequences of length 4 ([2,3,7,41] and [2,3,11,13]); three sequences of 
length 5 ([2,3,7,43,1805], [2,3,7,83,85] and [2,3,11,17,59]); 17 sequences 
of length 6; 27 sequences of length 7; and hundreds of sequences of length 8. 

So far, we know only 11 Giuga numbers (or proper Giuga sequences). They 
are 

3 factors: 

30 = 2 • 3 • 5 

4 factors: 

858 = 2 • 3 • 11 • 13 
1722 = 2 • 3 • 7 • 41 

5 factors: 

66198 = 2 • 3 • 11 • 17 • 59 

6 factors: 

2214408306 = 2 • 3 • 11 • 23 • 31 • 47057 
24423128562 = 2 • 3 • 7 • 43 • 3041 • 4447 

7 factors: 

432749205173838 = 2 • 3 • 7 • 59 • 163 • 1381 • 775807 
14737133470010574 = 2 • 3 • 7 • 71 • 103 • 67213 • 713863 
550843391309130318 = 2 • 3 • 7 • 71 • 103 • 61559 • 29133437 

8 factors: 
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244197000982499715087866346 = 

2 • 3 • 11 • 23 • 31 • 47137 • 28282147 • 3892535183 
554079914617070801288578559178 = 

2 • 3 • 11 • 23 • 31 • 47059 • 2259696349 • 110725121051 

For all of these examples, the ‘sum minus product’ value is 1; to reach any 
higher value, the sequence would have to have at least 59 factors. To find 
all Giuga sequences of a given length, one could check all sequences of this 
length whose elements are not too large (the sum over their reciprocals must 
be greater than 1). However, the number of these grows exponentially; even 
for length 7 there are too many to check them all. Fortunately, we have the 
following theorem which tells us how to find all Giuga sequences of length m 
with a given initial segment of length m — 2. 

Theorem 3. (a) Take an initial sequence of length m — 2, [m,..., n m _ 2 ]. 
Let 

P = n 1 • • • n m - 2 , S = 1 / 77,1 + • • • + l/n m _ 2 . 

Fix an integer v > S (this will be the sum minus product value). Take any 
integers a, b with a ■ b = P(P + S — v) and b > a. Let 

n m -1 := (P + a)/P(v - S), n m := (P + b)/P(y - S ). 

Then 

S + l/n m _i + l/n m - l/Pn m -±n m = v. 

The sequence [ni,... ,n m _i,n m ] is a Giuga sequence if and only if n m -\ is 
an integer. 

(b) Conversely, if \n \,... ,n m -i,n m ] is a Giuga sequence with sum minus 
product value v, and if we define 

a := n m -\P{v - S) - P, b := n m P{v - S) - P 

(with P and S the product and the sum of the first m — 2 terms) then a and b 
are integers and a ■ b = P(P + S — v). 

Proof. (a) First we have to check that with these definitions for n m _ 1 and 
n m we have in fact S + l/n m _i + 1 /n m — 1/Pn m _in m = v. These are 
straight-forward calculations: 

S + P(y - S)/(P + a) + P(v - S)/(P + b)~P 2 (v- S) 2 /P{P + a){P + b) = v 


10 



if and only if 


P/(P + a) + P/(P + b) - P{y - S)/(P + a)(P + b) = 1 

if and only if 

P/{P + a) + P/(P + b) - ( P 2 - ab)/(P + a)(P + b) = 1, 

which is true. 

This means that the completed sequence is a Giuga sequence if and only 
if both n m - 1 and n m are integers. It remains to be shown that n m _i is 
integer if and only if n m is. Because of the symmetry, it is enough to prove 
the implication in one direction. If 

P{y -S)\(P + a) 

then 

P(v - S) | {P + a){P + b) = 2 P 2 + {a + b)P - (v - S)P, 

so 

P(v - S) I (2 P 2 + {a + b)P) = P(P + a + P + b ), 

so 

P(v - S) I P{P + b) 

(since P(y — S ) divides P + a). 

The assertion follows if we show that gcd (P(v — S),P) = 1. Assume that 
there is a prime p with p \ P(v — S ) and p \ P, p \ rii, say. Since gcd(rij, rij) = 
1 for i not equal to j, p does not divide any of the other factors. Since p 
divides 

P{v - S) = vP - (n 2 • • • n m - 2 + ... + «!■•• n m _ 3 ), 
we can drop all terms on the right-hand side with a factor m to get 

p | Tl\ • • • Tli—l ' Tli+l ■ • • Tl m — 2, 

which is a contradiction. 

(b) It is clear that these a and b are integers. It remains to be checked that 
a ■ b = P(P + S — v). We have 

a-b = -P 2 (v - S)(n m -i+n m ) + P 2 n m -in m (v - S) 2 + P 2 

= P 2 (S - v)n m -in m ■ (l/n m _i + l/n m + S - v) + P 2 
= P 2 {S - v)n rn -in m ■ l/(n m -in m P) + P 2 
= P{P + S-v). © 
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We used this theorem to compute our Giuga sequences; this works well for 
sequences up to length 7, but for length 8 and higher there are still too many 
possible initial segments to check. 

All Giuga sequences we have found so far contain an even factor (it is usually 
the factor 2, but there are two sequences of length 9 which contain 4 instead). 
So far, we have not been able to find a sequence with only odd factors. 
Since any counterexample to Giuga’s conjecture would be an odd Giuga 
number, it would be of some interest to find at least one such sequence. Let 
n := n\ ■ ■ ■ n m . The Giuga equation l/ni + .. . + l/n m — l/n = v is equivalent 
to n/n\ + ... n/n m — 1 = nr; by considering this equation modulo 4, it is 
quite straightforward to show that if all factors n t are odd, then necessarily 
m — v = 1 mod 4. (In fact, assume that the first k factors are congruent 
to —1 modulo 4, and the other m — k factors are congruent to 1 modulo 4. 
Then the equation reduces to 

-1 = v(-l) k - j^-l)*" 1 - (m - fe)(-l) fc = (~l) k (v -m + 2k) mod 4, 

from which m — v = 1 mod 4 follows.) If we look for odd sequences with 
value v = 1, then we only have to check sequences of length m = 2 mod 4. 
The cases m = 2 and m = 6 can be ruled out because we need at least nine 
relatively prime odd integers for the sum of their reciprocals to exceed 1. 
Now, m = 10 can be ruled out computationally, with the use of Theorem 3. 
But this is where computational feasibility ends. For m = 14, there are just 
too many initial segments to check; another approach is needed here. 

We asserted in the introduction that there are infinitely many Giuga se¬ 
quences. As the following proposition tells us, it is possible to generate 
longer Giuga sequences out of shorter ones with certain properties. 


Theorem 4. Take a Giuga sequence of length rn, [m,..., n m _i, n m ], which 
satisfies 

1 T"m = Tl\ ' ' ' Tim— 1 1- (3) 


Let 

hm ■ ■ ■ n m —i + 1, h m+ i := n\ ■ ■ ■ n m —\fi m — 1. 

Then [n \,..., n m _i, h m , n m +i] is also a Giuga sequence with the same sum 
minus product value. 


Proof. Let P := n\... n m _i, S := 1 /ri\ + ... + Then n m = P — 1, 

h rn = P + 1 and h rn+ \ = P 2 + P — 1. Both sequences have the same sum 
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minus product value if and only if 


S + P~ 1 P(P- 1) S+ P+l + P 2 + P-l P{P+1){P 2 + P~ 1)’ 

the latter equation is true for all S and P. © 

Note that if the shorter Giuga sequence has property (3), then so has the 
longer one. Since the sequence [ 2 ,3, 5] also has this property, this proves 
that there are Giuga sequences of any length. (Note also that each of the 
sequences which occur in such a recursion will contain an even factor.) How¬ 
ever, the sequences arrived at by this recursion are not the only Giuga 
sequences there are. In each step from length m to m + 1, other Giuga 
sequences seem to pop out of thin air, some of them with property (3) (and 
thus leading to new recursions), some without. 

Explicitly, two infinite families are 

(a) m = 2 , n k = ni ■ ■ ■ n k -i + l for k = 2 ,..., nn - 1 , n m = n i • • • n m _i - 1 ; 

(b) n i = 2, ri 2 = 3, 713 = 11, 714 = 23, ns = 31, n k = n\ - ■ ■ n k ~ 1 + 1 for 
k = 6 ,..., to - 1 , n m = m ■ ■ ■ n m -1 - 1 - 

Recall that a Carmichael sequence is a finite increasing sequence of integers, 
[ai,..,, a m ], such that 

m \ 

[^[ (ij — 1 I for i = 1 ,..., m. 

=1 / 

At the end of the introduction we stated that there is a profusion of Carmichael 
sequences; in fact, there are infinitely many of them with any number of fac¬ 
tors. A trivial example would be the sequence [a,..., a] for any a € IN, but 
this is too cheap. We now conclude this section by giving a somewhat less 
trivial construction. We omit the computations here, since they are essen¬ 
tially simple but would enlarge this side remark unduly. 

There are infinitely many Carmichael sequences of length 3; in fact, the fol¬ 
lowing construction gives us all 3-factor Carmichael sequences. Take three 
integers bl, 62,63 € IN which are pairwise co-prime. Let c € JN be the solu¬ 
tion of the congruence (6162 + 6163 + 62 ^ 3 ) c = —(61 + 62 + 63 ) mod 6162 ^ 3 - 
(Such a solution exists and is unique modulo 616263 , because the integers 



13 



are pairwise co-prime. Equivalently, one can also solve the system of con¬ 
gruences & 2 & 3 C = —62 — b -3 mod 61 , 6163 c = —b\ — 63 mod 62 , b^c = —b\ — b -2 
mod 63 .) Then oq := cb\ + 1, <22 := cb 2 + 1, 03 := 063 + 1 is always 
a Carmichael sequence, and every 3-factor Carmichael sequence is of this 
form. For example, if we choose 61 = 1, 62 = 2 and 63 = 3, then we get 
Chernick’s observation (see [1]) that a\ := 6/s+l, <22 := 12 /c+l, 03 := 18/s+l 
is always a Carmichael sequence. 

The following two recursions produce Carmichael sequences of length m + 1 
and 772 + 2 out of the Carmichael sequence [<21, ..., a m ]. 


m 

(a) Let a m+ i := n aj. Then [cq,... ,a m ,a m+ 1 ] is a Carmichael sequence. 

3 1 

m 

(b) Let a m _|_i := n a j and a m+ 2 := c?(a m+ i — 1) + 1 where d is a divisor 

+ -1 

of a m+ i + 1. Then [ai,..., a m , a m+ i, a m+ 2 ] is a Carmichael sequence. 

In particular, if we start with an odd 3-factor Carmichael sequence, then we 
can extend it to arbitrary length by iterating the step (b) with d = 2. 


4 Open problems 

1. Giuga’s conjecture: Show that no integer exists which is both a Giuga 
number and a Carmichael number. 

More general: Show that no Giuga sequence can be a Carmichael se¬ 
quence. (This would imply the truth of Giuga’s conjecture.) 

2 . Does every Giuga sequence contain two factors n*, nj with nj \ (n^ —1)? 
If this were true, then Giuga’s conjecture is proved. 

3. Find a Giuga sequence which consists of odd factors (or odd primes) 
only, or prove that none exist. If there were none, then Giuga’s con¬ 
jecture would be proved. 

4. Are there infinitely many proper Giuga sequences? 

5. Is it true that, in the notation of Section 2, r m (L m ) = i m for m = 
5,..., 27692? If so, then this would prove that a counterexample to 
Giuga’s conjecture has at least 36069 digits. 
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6 . Find a fast way to generate all Giuga sequences of a given length. 

7. Are there Giuga sequences with a sum minus product value higher 
than 1 ? 

8 . Are there two distinct Giuga sequences whose elements have the same 
product? 

9. Can each integer be a factor in a Giuga sequence? If this were true 
then it would answer the previous two questions positively. In fact, 
take any integer n which is the product of a Giuga sequence, n = 
n\ ■ ■ ■ n m with 1/ni + ... + l/n m — 1/n = v. If we can find a sec¬ 
ond Giuga sequence which contains n as a factor, e.g., 1/n + 1/hi + 
... + 1/hk — 1/nhi • • • hfc = w, then we can combine the two of them 
and get the Giuga sequence l/n\ + ... + 1 /n m + l/h\ + .. .1/hk ~ 
1/ni • • • n m hi ■ ■ ■ hk = v + w. It has the same product as the previous 
one, but a higher sum minus product value. 

10. Agoh’s conjecture: Let Bk denote the kth Bernoulli number. Then 
nB n _ i = —1 mod n if and only if n is a prime? Note that the 
denominator of the number nB n -i can be greater than 1 , but since the 
denominator of any Bernoulli number is squarefree, the denominator 
of nB n _ 1 is invertible modulo n. As Takashi Agoh (Science University 
of Tokyo) has informed us, this recent conjecture of his is equivalent 
to Giuga’s conjecture: Every counterexample to Giuga’s conjecture is 
also a counterexample to Agoh’s conjecture and vice versa. This can 
be seen from the well-known formula 

Sn-i = n l ~ 1 B ri ^ l 

after some analysis involving von Staudt-Clausen’s theorem: The de¬ 
nominator of B 2k is given by Jj P- (See [4], pp. 91-93; this also 

p prime 
(p-l)|2fc 

implies that the denominator of B 2 k is squarefree.) 

Incidentally, it is possible to use a similar argument to characterize 
Giuga numbers in the following way: n is a Giuga number if and only 
if nB^ n ) = -1 mod n. 

Finally, we would like to thank Hugh Edgar for originally making us aware 
of Giuga’s conjecture and challenging us to extend what was known compu¬ 
tationally. 
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